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Abstract — The results of photometric and polarimetric observations carried out for some bright atmosphere - 
less bodies of the Solar system near the zero phase angle reveal the simultaneous existence of two spectacular 
optical phenomena, the so-called brightness and polarization opposition effects. In a number of studies, these 
phenomena were explained by the influence of coherent backscattering. However, in general, the interference 
concept of coherent backscattering can be used only in the case where the particles are in the far-field zones 
of each other, i.e., when the scattering medium is rather rarefied. Because of this, it is important to prove rig- 
orously and to demonstrate that the coherent backscattering effect may also exist in densely packed scattering 
media like regolith surface layers of celestial bodies. From the results of the computer modeling performed 
with the use of numerically exact solutions of the macroscopic Maxwell equations for discrete random media 
with different packing densities of particles, we studied the origin and evolution of all the opposition phenom- 
ena predicted by the coherent backscattering theory for low-packing- density media. It has been shown that 
the predictions of this theory remain valid for rather high packing densities of particles that are typical, in par- 
ticular, of regolith surfaces of the Solar system bodies. The results allow us to conclude that both opposition 
effects observed simultaneously in some high-albedo atmosphereless bodies of the Solar system are caused 
precisely by coherent backscattering of solar light in the regolith layers composed of microscopic particles. 

DOI: 10. 1 134/S0038094613060014 


INTRODUCTION 

The results of photometric and polarimetric obser- 
vations carried out near the zero phase angle show that 
two unique optical phenomena, the so-called bright- 
ness and polarization opposition effects, are simulta- 
neously present in some bright atmosphereless bodies 
of the Solar system. Figure 1 taken from the paper by 
Mishchenko et al. (2006a) displays the observational 
data acquired in the visible spectral range for the high- 
albedo asteroids 44 Nysa and 64 Angelina, the Jovian 
satellite Europa, and Saturn’s rings (Franklin and 
Cook, 1965; Zellner and Gradie, 1976; Johnson et al., 
1980; Harris et al., 1989; Thompson and Lockwood, 
1992; Dollfus, 1996; Rosenbush et al., 2005; Rosen- 
bush and Kiselev, 2005). It is seen that each of the pho- 
tometric phase curves in the left part of the plot has a 
very narrow nonlinear peak centered at zero phase 
angle and superimposed on a linear background. This 
nonlinear peak of brightness is called the brightness 
opposition effect (BOE). In the right part of the plot, 
each of the phase dependencies of the degree of linear 
polarization demonstrates a narrow local minimum of 
negative polarization centered at a phase angle 
approximately equal to a half-width of the corre- 
sponding brightness peak. Each minimum is super- 
posed on a broad, nearly parabolic branch of negative 
polarization that is characteristic of the majority of 
atmosphereless bodies of the Solar system. Such a nar- 


row asymmetric minimum of negative polarization is 
called the polarization opposition effect (POE) 
(Mishchenko, 1993). 

Shkuratov (1989) and Muinonen (1990) were the 
first who hypothesized that two optical effects 
observed in a wide class of atmosphereless bodies near 
opposition, namely BOE and the parabolic branch of 
negative polarization, are caused by coherent back- 
scattering in discrete random media. However, it is 
known that the parabolic branch of negative polariza- 
tion can be induced by both single scattering by non- 
spherical grains of the surface regolith layer and near- 
field effects (Petrova and Tishkovets, 2011). In the 
papers by Mishchenko (1993) and Mishchenko et al. 
(2000), it was concluded that coherent backscattering 
causes an extremely narrow asymmetric peak of negative 
polarization that appears if the incident radiation is non- 
polarized. It was supposed that the BOE should be 
accompanied by the POE, if a substantial portion of the 
surface is covered by particles smaller than the wave- 
length, and the angle of the minimal negative polariza- 
tion should be comparable to the BOE half-width. 

In the papers by Mishchenko (1992a; 1992b), the 
rigorous relations allowed all the characteristics of the 
coherent backscattering effect to be calculated for the 
precise backscattering direction from the exact solu- 
tion of the vector radiative transfer equation for the 
low-packing-density media composed of particles of 
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Fig. 1 . The values of the normalized (to the minimal phase angle) intensity and the degree of linear polarization of light measured 
in the asteroids 44 Nysa and 64 Angelina, the satellite Europa, and Saturn’s rings (Mishchenko et al., 2006a). 


arbitrary shape and size. With the use of these rela- 
tions, it was shown that the coherent backscattering 
effect in the regolith layers containing submicron par- 
ticles with a refractive index of about 1.31 (for the sat- 
ellite Europa and Saturn’s rings) and 1.65 (for the 
high-albedo asteroids 44 Nysa and 64 Angelina) can 
be one of the causes of the observed brightness opposi- 
tion effects (Mishchenko and Dlugach, 1992; 1993; 
Dlugach and Mishchenko, 1999). Moreover, this the- 
ory was successfully used (Mishchenko and Dlugach, 
2008) in the interpretation of the radar polarization 
observations of Saturn’s rings, namely, the results of 


measurements of the circular polarization ratios 
(Ostro et al., 1980; Nicholson et al., 2005). 

However, the interference concept of the coherent 
backscattering effect is strictly applicable only when 
the particles are in the far scattering zone relative to 
each other, i.e., when the scattering medium is suffi- 
ciently sparse (the packing density of particles is neg- 
ligible) (Barabanenkov et al., 1991; Mishchenko et al., 
2006b). Therefore, one faces a question of whether the 
coherent backscattering effect can also take place in 
the case of densely packed media like regolith layers of 
the surfaces of celestial bodies. A definitive answer to 
this question can only be obtained from numerically 
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Fig. 2. Schematic explanation of the physical nature of the brightness and polarization opposition effects. 


exact solutions of the macroscopic Maxwell equations 
for the media composed of varying numbers of ran- 
domly positioned particles. This approach allows, 
first, the uncertainties connected with the use of dif- 
ferent approximations to be eliminated and, second, 
the influence of different parameters of the medium to 
be immediately controlled by individual varying of 
their values. Thus, one may study the evolution of all 
effects caused by coherent backscattering as the pack- 
ing density of particles grows from zero to the values 
typical of the regolith surface layers, and the results 
obtained allow of a clear answer on the physical nature 
of the opposition effects observed. 

THEORY 

The effect of coherent backscattering (or weak 
localization) of electromagnetic waves in discrete ran- 
dom media predicted by Watson ( 1 969) in the study of 
multiple scattering of electromagnetic waves in a rar- 
efied plasma still remains the subject of active theoret- 
ical and laboratory researches (see, e.g., Akkermans 
et al., 1988; van der Mark et al., 1988; Barabenkov et al., 
1991; Kuz’min and Romanov, 1996; Mishchenko et al., 
2006b). To explain the interference nature of the 
coherent backscattering effect, let us consider a dis- 
crete medium composed of randomly positioned scat- 
tering particles illuminated by a plane wave incident in 
the direction specified by the unit vector n m (Fig. 2a). 
Let the direction of the scattered wave be determined 
by the unit vector n obs , and the observer be very far 
from the medium. Consider two reverse (conjugate) 
scattering paths involving the same configuration of 
Af particles. The waves scattered through these paths 
will interfere, and the interference is constructive or 
destructive depending on the phase difference between 
the paths. If the observer is sufficiently far from the 
exact backscattering direction (-n m ), the average 
effect of interference of the conjugate scattered waves 
going through a group of particles in opposite direc- 
tions is zero, owing to the randomness of particle posi- 
tions. Consequently, the observer measures some 


incoherent (diffuse) intensity. However, at a = 0° 
(ii obs = -n m ), the phase difference between the conju- 
gate paths involving any configuration of particles is 
identically equal to zero, i.e., the coherence com- 
pletely remains, and the interference is always con- 
structive. It is this effect that induces the opposition 
peak in the intensity of the scattered radiation. 

To explain the nature of the polarization opposition 
effect, let us consider particles 1—4 lying in a plane 
normal to the illumination direction (Fig. 4b) and 
assume that their sizes are much smaller than the 
wavelength (Shkuratov, 1989; Muinonen, 1990). Parti- 
cles 1 and 2 lie in the scattering plane, while particles 3 and 
4 are in the perpendicular plane. If the incident radia- 
tion is unpolarized, then it can be represented as a 
mixture of two linearly polarized beams of light with 
mutually perpendicular directions of the polarization 
planes. Light scattered along conjugate paths involv- 
ing particles 3 and 4 is negatively polarized, while scat- 
tering along the paths involving particles 1 and 2 yields 
positive polarization. In the first case, the phase differ- 
ence between the conjugate paths is always zero, while 
in the second case it is zero only at a = 0° and oscil- 
lates rapidly with increasing a. Therefore, on average, 
coherent backscattering enhances the contribution of 
negative polarization over a wider range of phase 
angles. The result is the negative polarization mini- 
mum at a small phase angle (POE) comparable to the 
angular width of the coherent intensity peak (BOE). 
The fact that only certain particle configurations con- 
tribute to the POE often makes the latter less pro- 
nounced than the BOE. 

Apparently, the POE was observed for the first time 
in laboratory by Lyot (1929), when he measured the 
degree of polarization of a powder-like surface com- 
posed of microscopic MgO particles. The author 
called the measured narrow minimum of negative 
polarization mysterious and attributed its existence to 
extremely small sizes of magnesium oxide particles. In 
the study by Shkuratov et al. (2002), these results were 
successfully reproduced and supplemented by photo- 
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metric measurements that showed the simultaneous 
presence of a narrow BOE. 

The interference explanation of the coherent back- 
scattering effect relies on the notion of phase of an 
electromagnetic wave. In this connection, it is impor- 
tant to stress that the notion of phase itself is applicable 
only to transverse electromagnetic waves, such as a 
plane or spherical wave. Due to this, the interference 
explanation of the coherent backscattering effect indi- 
rectly relies on the assumption that each of the parti- 
cles in any sequence of particles (Fig. 2a) is in the far- 
field scattering zones of both the previous particle and 
the next one. Therefore, in order to justify the involve- 
ment of the concept of coherent backscattering for 
explaining the observational data acquired for densely 
packed media, the range of applicability of the coher- 
ent backscattering theory should be clearly deter- 
mined. This can be done on the basis of the results of 
the numerically exact solution of the macroscopic 
Maxwell equations for scattering media with a varied 
packing density of particles. Such computations have 
become possible only recently due to the development 
of effective numerical methods and the increase of the 
computer power (especially, due to the ever-increasing 
power of computers (Mackowski and Mishchenko, 
2011; and references therein)). 


THE TECHNIQUE OF NUMERICALLY EXACT 
MODELING OF LIGHT SCATTERING 
BY A MEDIUM WITH AN ARBITRARY 
PACKING DENSITY OF PARTICLES 

To examine with numerically exact calculations 
whether the use of the coherent backscattering theory 
is valid for densely packed media, an adequate model 
of the scattering medium should, first of all, be chosen. 
As such a model, let us consider an imagined spherical 
volume of radius R filled with TV randomly positioned 
nonoverlapping spherical particles of radius r (see Fig. 3). 
The size parameter of the scattering volume and each 
of the particles is kR and kr , respectively, where k is the 
wave number in vacuum. To model the statistically ran- 
dom particle positions within the spherical volume con- 
sistent with the fundamental assumption of ergodicity 
(i.e., for a specified realization of a random scattering 
process, the average over time is equivalent to the average 
over an ensemble (Mishchenko et al., 2006b)), we follow 
the approach pioneered by Mishchenko et al. (2007). 
Specifically, we use one realization of an TV-particle 
group generated randomly according to the algorithm 
by Mackowski (2006) and then average all scattering 
characteristics over the uniform orientation distribu- 
tion of this configuration with respect to the laboratory 
coordinate system. The result is an infinite continuous 
set of random realizations of the TV-particle group, 
which enables one to use the efficient orientation aver- 
aging technique developed in the framework of the 
superposition F-matrix method (Mackowski and 
Mishchenko, 1996; Mishchenko et al., 2002). 



Fig. 3. Schematic explanation of the light scattering pro- 
cess in a macroscopic spherical volume filled with Anon- 
overlapping randomly positioned spherical particles. 


It is obvious that the adopted model cannot repro- 
duce the infinite diversity of the particulate media 
morphologies encountered in natural conditions and 
in the laboratory. However, it proves to be sufficiently 
representative to permit a clear pattern of the effects of 
packing density on light scattering by a discrete ran- 
dom medium. 

Assume that the statistically random particulate 
volume is illuminated by a parallel quasi-monochro- 
matic beam of light propagating in the direction of the 
unit vector n m and the observer is located in the far- 
field zone of the entire scattering volume in the direc- 
tion n obs (Fig. 3). Since the scattering properties of the 
volume are averaged over all orientations of the /V- par- 
ticle random configuration, they depend only on the 
scattering angle 0 (or the phase angle a = n — 0) pro- 
vided that the scattering plane is used for defining the 
Stokes parameters of the incident and scattered light. 

In the far-field zone, the Stokes parameters for 
the incident and scattered radiation are related as 
(Mishchenko et al., 2006b) 
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Fig. 4. The scattering characteristics calculated for spherical volumes randomly filled with spherical particles. 
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The scattering matrix F(a) is a special case of the 
Mueller transformation matrix (Hovenier et al. , 2004) . 
Strictly speaking, the elements of the matrix F(a) 
denoted in Eq. (1) by zeros vanish only if the primor- 
dial TV- sphere configuration had a plane of symmetry 
Our calculations showed, however, that in all cases 
considered below the magnitudes of these matrix ele- 
ments are much smaller than those of the other ele- 
ments in absolute value. In other words, the zeros 
denote the matrix elements that are negligibly small 
(in an absolute sense) relative to the other elements at 
the same scattering angles. 

From the scattering matrix elements, the conven- 
tional optical characteristics, corresponding to differ- 
ent types of polarization state of the incident radia- 
tion, can be calculated. Specifically, if the incident 
light is unpolarized, the first element a x ( a) (the phase 
function) determines the angular distribution of the 
intensity of the scattered radiation in the far field, 
while the ratio — Z? 1 (a)/a 1 (a) gives the corresponding 
degree of linear polarization. If the incident radiation 
is polarized linearly in the scattering plane, the linear 
polarization ratio is defined as 

= fli(a)-g?(q) (2) 

L a) + 2b x (a) + a 2 ( a) 

If the incident radiation is polarized circularly in the 
counterclockwise direction, when looking in the 
direction of propagation, the circular polarization 
ratio is 

= gi(a) + a 4 (a) ^ 

a { (a) - a 4 (a) 

It is worth stressing that the degree of linear polar- 
ization of light singly scattered by one particle (for 
unpolarized incident light) may have a negative branch 
in the range of phase angles from 0° to 20°. In such 
cases, a direct computer modeling of the POE can be 
quite problematic, since it is impossible to distinguish 
between the singly and multiply scattered contribu- 
tions in a direct numerical solution of the Maxwell 
equations (Mishchenko et al., 2007; Petrova et al., 
2007). In order to eliminate this uncertainty, the fol- 
lowing technique was suggested (Mishchenko et al., 
2009a). For each value of the considered refractive 
index, we select a size parameter value kr yielding a 
single-scattering contribution with near-zero values in 
a wide interval at small phase angles. This allows the 
multiple -scattering contribution to be unambiguously 
distinguished and qualitatively estimated. Note that 
this approach made it possible to demonstrate for the 
first time, with the use of numerically exact computer 
solutions of the Maxwell equations, the whole set of 
optical effects predicted by the low-density theory of 
coherent backscattering (Mishchenko et al., 2009a; 
2009b; Dlugach et al., 201 1). 

Different optical characteristics of the scattering 
volumes containing both nonabsorbing and absorbing 


particles were calculated with the superposition F- matrix 
method (Mackowski and Mishchenko, 1996; 2011). The 
cases of polydisperse scattering media composed of 
particles of different sizes or different chemical com- 
position were also considered. Here, we will only 
briefly review the results obtained for the scattering 
volumes with the size parameter kR ranging from 20 to 
60. The volumes are filled with nonabsorbing spherical 
particles; several values of the real refractive index 
were considered: m = m r = 1.194 (latex in water), 1.31 
(water droplets or ice in the visible spectral range) ,1.4, 
1.5, and 1.6 (different minerals). Depending on the 
value of kR , the number of particles TV was varied; and 
TV max = 900 for kR = 30. 

Figure 4a illustrates the results of our calculations 
for the phase range 0° < a < 30°. They were performed 
for m = 1 .3 1 , kR = 40, kr = 2, and the packing density p, 
changing from 1.5% (N = 100) to 11.7% (N = 800). 
For comparison, we also show the data corresponding 
to the scattering by a single particle. It is seen that the 
results of a numerically exact solution of the Maxwell 
equations completely agree with the qualitative pre- 
dictions of the low-density theory of coherent back- 
scattering, namely: 

(1) The normalized intensities of the scattered 
radiation a^aj/a^O 0 ) exhibit the backscattering peaks 
rapidly developing with TV, and the angular width of 
these peaks remains practically the same. 

(2) In the phase-angle range from 0° < a < 30°, the 
degree of linear polarization -b^ofy/a^a) is equal to 
zero for TV = 1, but rapidly develops a pronounced 
minimum with growing TV caused by the increasing 
contribution from multiple -scattering; the phase angle 
of the polarization minimum a min is virtually indepen- 
dent of TV and is comparable to the angular width of the 
peak of the normalized intensity a x { dj/a^O 0 ) . Further- 
more, the angular shape of the branch of negative 
polarization is asymmetric, with being signifi- 
cantly closer to zero than to the inversion phase angle. 

In order to analyze the dependence of the optical 
characteristics on the size parameter of the scattering 
volume kR, the results of calculations for m = 1.31, 
N= 800, and kR = 30, 40, and 60 are shown in Fig. 4b. 
It is seen that the angular width of all peaks and the 
angle of the polarization minimum decrease with 
increasing kR approximately as 1/kR. 

In Fig. 4c, the results of calculations for the con- 
stant value of the size parameter kR = 40 and the 
monomer number TV = 500, and for the refractive 
indexes m = 1.31, 1.4, 1.5, and 1.6 are presented. One 
can easily see similar opposition effects for all refrac- 
tive index values, which illustrates the interference 
nature of the coherent backscattering effect. 

From Figs. 4a— 4c, we may conclude that all of the 
peculiarities in the behavior of the peaks of normalized 
intensity and polarization minimum confirm their 
interference origin caused by multiple scattering. Spe- 
cifically, the peaks of intensity of the scattered radia- 
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tion and the minima of negative polarization are 
absent in the case of a single particle, but appear for 
the scattering volumes with N > 1. The angular width 
of the normalized intensity peak is comparable to the 
phase angle a min , and both of them are independent of 
TV and m and are inversely proportional to the size 
parameter of a scattering volume kR. Furthermore, 
the polarization minima become increasingly asym- 
metric with increasing kR. The opposition peaks of 
intensity appear and develop with increasing TV due to 
the influence of the multiple -scattering contribution, 
and their angular widths are inversely proportional to 
kR , which once again confirms their interference ori- 
gin. For a fixed kR , the angular widths of the peaks are 
approximately the same and are independent of the 
number of particles and their refractive index, though 
the peak amplitudes increases with increasing refrac- 
tive index. 

The linear p L and circular p c polarization ratios 
calculated with Eqs. (2) and (3) can serve as an even 
better illustration of the origin and evolution of the 
coherent backscattering effect, since they coincide 
with the abscissa axis for a single particle. The back- 
scattering peaks appear and grow with increasing TV 
under the influence of the multiple -scattering contri- 
bution, and their angular widths are inversely propor- 
tional to kR , which is indicative of their interference 
origin. For a fixed kR , the angular width is practically 
independent of the number of particles and their 
refractive index, though the peak height increases with 
increasing refractive index. 

The effect of increasing the number of particles in 
the scattering volume can be expected to be twofold. 
Initially, it simulates multiple scattering and thus 
enhances the classical manifestations of diffuse and 
coherent scattering of light. Eventually, however, the 
increased number of particles can cause changes in the 
scattering characteristics of the whole volume not 
implied by the classical theories of radiative transfer and 
coherent backscattering. One should, therefore, expect 
that the analysis of the numerically exact T-matrix results 
with the use of these theories becomes inadequate, when 
the packing density of particles increases to rather large 
values. 

In order to analyze the effect of increasing the 
number of particles TV on the scattering properties of 
the volume in the backscattering direction, we con- 
sider the results of calculations carried out for kR = 20, 
m = 1.31, kr= 2, and TV varied from 50 (p ~ 7%) to 345 
(p ~ 47%) (Fig. 4d). We see that, indeed, (as the calcu- 
lations showed, the effect appears starting from TV = 
300 (N = 300 (p ~ 41%)) the curves develop interfer- 
ence ripples typical of a single spherical particle with a 
radius substantially larger than the wavelength (calcu- 
lations performed for kR = 31 show a similar result 
(Mackowski and Mishchenko, 2011)). Obviously, 
such a behavior of all of the considered optical charac- 
teristics is not predicted by the coherent backscatter- 


ing theory valid for low-density media. Nevertheless, 
the results of our computations do demonstrate that the 
predictions of this theory can survive relatively high pack- 
ing densities of particles typical, in particular, of the 
regolith surfaces of Solar system bodies. It is worth noting 
that the recent results (Muinonen et al., 2012) 
obtained with the Monte-Carlo technique for calcula- 
tions of the diffuse and coherent components accord- 
ing to the microphysical theory of radiative transfer 
and coherent backscattering (Mishchenko et al., 
2006b) confirm this statement. 

Our results of the numerically exact modeling per- 
formed for the media with different packing densities 
pose a question as to why the manifestations of the 
coherent backscattering effect are so remarkably 
immune to the increase of the packing density of par- 
ticles. This can be caused by the fact that, even in rel- 
atively densely packed media, the partial wavelets scat- 
tered by sequences of widely separated particles still 
provide a sufficient contribution to the total scattered 
radiation and, thereby, make the classical coherent 
backscattering features quite pronounced. 

CONCLUDING REMARKS 

It must be emphasized once again that the results of 
computations presented above were obtained for a rel- 
atively simple model of a scattering medium. More- 
over, one should have in mind that the disk-integrated 
observations of the celestial bodies having the hetero- 
geneous surfaces yield cumulative information on the 
contribution of morphologically different surface 
areas with different values of albedo. Some of them 
may induce the BOE and POE with varying angular 
width and amplitude, while the others may cause only 
the BOE changing over the surface. These factors 
make the prospects of unambiguous interpretation of 
such observational data, especially with the direct 
solution of the Maxwell equations in terms of the 
refractive index, size, and packing density of particles, 
to be extremely problematic. Nevertheless, the possi- 
bility of the presence of the coherent backscattering 
effect in rather densely packed media confirmed by 
our modeling allows the opposition effects observed in 
a number of high-albedo atmosphereless bodies of the 
Solar system (Fig. 1) to be reliably explained. 

The calculated values of the opposition peak 
amplitude in intensity of radiation and the opposition 
minimum depth in polarization (Figs. 4a, 4c) are 
larger than the corresponding values measured in the 
high-albedo asteroids and planetary satellites (Rosen- 
bush and Mishchenko, 2011). This difference would 
be expected, since it is highly improbable that the 
entire surface of these celestial bodies is covered by 
identical microscopic grains inducing the spatially 
constant BOE and POE. The measured angular width 
of the opposition effects turns out to be substantially 
smaller than the calculated one, which can be 
explained by the fact that the interference base for a 
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finite scattering volume is determined by its size 
parameter kR , while for the optically thick nonabsorb- 
ing or weakly absorbing regolith layer it is controlled by 
the product kX tr , where A, tr is the transport free path 
(Barabanenkov et al., 1991). The actual values of kX ir 
may substantially exceed the kR values used in our 
computations, thereby resulting in much narrower 
opposition effects observed. Note that the Monte- 
Carlo calculations (Muinonen et al., 2012) showed 
that the substantial increase of the size of a scattering 
volume really causes the considerable decrease of the 
angular width of the opposition peaks. 

In conclusion, we stress that the photometric and 
polarimetric observations of the asteroids 44 Nysa and 
64 Angelina, the satellite Europa, and Saturn’s rings 
(Fig. 1) are unique in that they reveal simultaneous 
existence of the brightness and polarization opposition 
effects of nearly equal angular widths and angular pro- 
files consistent with the exact solutions of the Maxwell 
equations (Fig. 4). Up to now, no other theory directly 
based on the Maxwell equations, except the coherent 
backscattering theory, has predicted both these effects 
simultaneously. Therefore, our numerically exact 
results allow of the conclusion that both opposition 
effects are caused by coherent backscattering of solar 
light in the regolith layers composed of microscopic 
ice (for the satellite Europa and Saturn’s rings) and sil- 
icate (for the asteroids 44 Nysa and 64 Angelina) par- 
ticles. 
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